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ABSTRACT 


A  high-current,  high-power  X-ray  Free  Electron  Laser  (FEL)  using  the  Stanford 
linear  accelerator  (SLAG)  as  an  electron  beam  source  has  been  proposed.  Such  a 
system  will  provide  the  scientific  community  with  its  first  realistic  X-ray  laser.  One 
dramatic  use  a  hard  X-ray  laser  light  at  the  proposed  wavelength  of  1.5  Angstrom 
would  be  to  image  DNA  base-pairs  [7]. 

With  an  accelerator  already  in  place,  the  major  cost  of  developing  the  SLAG 
system  would  be  in  the  manufacturing  of  the  undulator.  The  cost  in  building  an 
undulator  is  proportional  to  its  length.  This  thesis  uses  numerical  simulation  to 
evaluate  the  effectiveness  of  an  FEL  klystron  in  order  to  decrease  the  length  of  the 
undulator  with  a  dispersive  section. 

Simulations  show  that  the  quality  of  the  electron  beam  injected  into  the  undulator 
can  greatly  effect  the  ability  of  an  FEL  to  produce  coherent  light,  and  is  a  factor  in 
determining  the  length  of  an  undulator  needed  for  saturation.  In  the  high-gain  FEL 
(such  as  proposed  at  SLAG),  a  dispersive  section  is  sensitive  to  energy  spread  and 
may  adversely  affect  the  gain  of  the  system.  For  moderate  to  low  energy  spread,  it  is 
found  that  a  klystron  is  useful  in  shortening  the  undulator  length. 

By  examining  how  the  strength  of  a  klystron  dispersive  section  affects  the  gain  of 
a  high-current  system,  we  determine  a  critical  strength  above  which  the  klystron  can 
be  detrimental.  A  dispersive  section  that  is  either  too  strong  or  not  strong  enough  will 
result  in  less  than  optimum  gain,  or  an  increase  in  the  required  length  of  an  undulator. 

Single-mode,  phase-space  simulations  are  used  to  investigate  the  effect  on 
electron  bunching  and  the  onset  of  saturation.  Longitudinal  multimode  simulations 
show  the  resulting  coherence  development. 
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I.  INTRODUCTION 


With  the  present  world  turmoil,  the  political  and  economic  structure  of  current 
nations  has  made  it  easier  for  the  Third  World  nations  or  even  small  factions  of 
terrorist  organizations  to  obtain  state-of-the-art  weapons.  The  advances  in  technology 
as  well  as  hardening  of  modern  day  missiles  has  driven  present  day  ship  self-defense 
systems  obsolete.  The  time  required  to  track,  engage,  and  ensure  destruction  of 
incoming  weapons  has  decreased  to  seconds.  Present  day  defense  systems  can  no 
longer  protect  ships  adequately. 

The  present  day  ship  self-defense  system  is  a  layered  system  that  relies  on  long 
and  medium  range  missiles  for  long  distance  engagements.  A  ship’s  gun  and  short 
range  missiles  such  as  RAM  and  SEA  SPARROW  are  for  missiles  that  get  through  the 
long  range  defenses  and  are  at  a  medium  range.  The  final  and  most  crucial  system  in 
the  layered  defense  is  the  Phalanx  Close-ln-Weapon  System  (CIWS),  for  close  in 
missiles  threats.  All  of  these  systems  except  for  CIWS.  are  not  always  available  to 
every  platform.  Even  if  the  outer  defense  systems  are  available,  the  detection  and 
engagement  of  these  incoming  missiles  at  long  ranges  have  become  difficult  with  the 
design  of  sea  skimming  and  low,  slow  flying  missiles.  The  present  self-defense 
systems  must  also  overcome  a  large  scale  attack  with  multiple  weapons  being  fired 
upon  them  simultaneously.  A  ship’s  ability  to  protect  itself  against  such  an  attack  is  at 
present  unrealistic. 

CIWS  has  a  reliable  detection  range  of  about  5000  m  and  an  effective 
engagement  range  of  less  than  2000  m.  A  typical  missile  requires  numerous  hits  from 
CIWS  before  it  is  destroyed.  With  the  missiles  traveling  at  speeds  up  to  and  greater 
than  mach  3,  the  range  that  CIWS  actually  kills  the  missile  is  much  closer  than  the 
2000  m  engagement  range.  Typically,  the  actual  kill  range  is  around  600  m.  If  CIWS 
is  able  to  destroy  the  missile,  fragments  from  the  missile  are  still  large  enough  and 
traveling  with  enough  momentum  to  reach  the  ship  and  create  significant  damage.  If 
the  incoming  missile  Is  designed  to  maneuver  once  it  is  in  close  range  of  its  target,  a 
CIWS  would  have  to  slew  its  gun  barrels  causing  the  rounds  to  be  off  target  until  the 
missile  stabilizes  and  the  CIWS  is  locked  on  again.  The  time  that  a  CIWS  round 
travels  to  reach  the  incoming  missile  also  adds  to  the  decreasing  range  problem.  It  is 
evident  the  technology  of  missiles  has  outgrown  ships’  self-defense  systems. 
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There  is  need  for  a  speed-of-light  weapon  which  can  reach  farther  than  CIWS 
and  kill  an  incoming  missile  faster,  so  as  to  allow  for  multiple  threat  engagements.  To 
kill  an  incoming  missile  with  a  laser,  the  power  required  on  the  surface  of  the  missile  is 
dictated  by  the  Fluence,  F  =  PAtiA  [1].  Fluence  (in  Jlcm^),  is  the  common  measure 
of  "hardness"  of  a  missile.  It  describes  the  amount  of  energy  absorbed  by  the 
missile’s  skin.  A  moderately  hardened  missile  might  require  a  fluence  of  10  kJlcm^ 
[1].  A  spot  size  A  =  5.0  cm^  is  achieved  at  a  range  R  =  5.0  km,  wavelength 
A,  =  3.8  microm ,  with  a  laser  aperture  of  D  ==  1 .5  meters,  where 

P  is  the  power  received  at  the  missile  and  Af  is  the  dwell  time  required  to  kill  the 
missile.  For  a  0.1  second  dwell  time,  the  total  power  to  destroy  the  missile  is 
approximately  P  =  500  kW.  The  ship’s  weapon  must  emit  greater  power  than  this  to 
overcome  aerosols  at  sea  level  over  the  distance  it  travels.  The  power  transmitted 
through  the  atmosphere  decreases  exponentially  as  P(z)  =  Po®^>  where 
a  =  5  X  10"^  km~''  is  the  extinction  coefficient  for  a  wavelength  of  A,  =  3.8  microns,  Pq 
is  the  initial  power  at  the  source,  and  z  is  the  range  to  the  target.  At  a  r^ge  of  20  km 
the  power  required  at  the  source  is  about  1 .4  MW  or  greater.  This  amount  of  power 
has  been  achieved  with  chemical  lasers.  The  minimal  amount  of  time,  ranging  from 
microseconds  to  seconds,  that  the  missile  will  travel  while  the  laser  is  locked  on, 
greatly  decreases  the  closing  distance  the  missile  travels  to  only  meters.  This  is  a 
significant  difference  compared  to  the  distance  the  missile  closes  in  on  its  target 
before  a  CIWS  obtains  a  kill,  and  it  is  performed  at  a  range  that  is  much  greater  than 
the  engagement  range  of  CIWS.  The  laser’s  speed  to  target  also  helps  to  reduce  the 
closure  distance.  The  range  at  which  the  missile  is  hit,  even  if  not  destroyed  but 
damaged,  is  sufficient  to  prevent  the  incoming  weapon  from  ever  reaching  it’s  target. 

A  Free  Electron  Laser  (FEL)  has  many  advantages  over  conventional  lasers  in 
producing  high  power  outputs.  Conventional  lasers  of  sufficient  power  use  dangerous 
chemicals  that  are  difficult  to  store  and  are  partially  depleted  with  every  firing  of  the 
laser.  The  depletion  of  the  source  for  lasing  then  limits  the  number  of  times  you  can 
fire,  and  the  number  of  incoming  missiles  that  can  be  engaged.  An  FEL’s  power 
source  is  the  electrical  power  from  the  ship,  which  is  degraded  only  if  the  ship  is 
already  badly  damaged.  The  FEL  is  tunable  over  a  wide  band  of  frequencies.  That 
allows  you  to  compensate  for  atmospheric  conditions  and  concentrate  the  light  spot  on 
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the  inbound  missile. 

The  Free  Electron  Laser  has  other  significant  uses  that  can  be  exploited  by  the 
military,  other  than  ship’s  self-defense.  There  are  many  uses  in  medicine  such  as 
DNA  mapping  and  X-rays.  Target  tracking  and  identifying  can  be  enhanced  with  a 
tunable  laser.  Ideally,  the  burden  of  research  and  development  can  be  shared  by 
civilian  and  military  ensuring  continuous  upgrading  and  off-the-shelf  products. 


II.  BACKGROUND 


A.  FEL  THEORY 

There  are  two  major  requirements  for  a  FEL  to  produce  emission;  a  relativistic 
electron  beam  supplied  by  an  accelerator,  and  an  undulator  with  a  spatial  periodic 
magnetic  field  [2].  The  FEL  creates  coherent,  monochromatic  light  by  accelerating 
relativistic  electrons  in  the  transverse  direction.  The  periodic  magnetic  field  from  the 
undulator  exerts  a  force  on  the  relativistic  electron  beam  as  it  passes  through  the 
undulator  causing  the  electrons  to  "wiggle"  in  the  transverse  direction.  The  electrons 
are  accelerated  from  side  to  side  causing  them  to  spontaneously  emit  light  in  a 
forward  cone  along  the  electron  beam  axis.  To  collect  a  fraction  of  the  spontaneously 
radiated  electromagnetic  field,  mirrors  are  placed  on  the  beam  axis  beyond  both  ends 
of  the  undulator  creating  an  optical  cavity.  The  coupling  created  by  passing  the 
electron  beam  through  the  stored  electromagnetic  field  causes  stimulated  emission. 
As  a  result,  the  optical  field  grows  and  coherent  light  is  produced  from  the  FEL  [2]. 

The  electrons  are  traveling  relativistically  with  energy  ymc^  prior  to  entering  the 
undulator,  where  m  is  the  electron  rest  mass  and  c  is  the  velocity  of  light.  The 
Lorentz  factor  of  the  electron  is  y=  where  p  =  vie  and  v  is  the  magnitude  of 

the  electron  velocity.  Upon  entering  the  undulator,  a  force  from  the  magnetic  field  acts 
upon  the  electrons  causing  the  electrons  to  move  back  and  forth  in  the  transverse 
direction. 

The  optical  field  properties  can  be  controlled  by  controlling  the  parameters  of  the 
undulator.  An  undulator  can  be  of  different  designs  using  either  permanent  magnets 
or  electromagnets  with  either  linearly  or  circularly  polarized  magnetic  fields  producing 
linearly  or  circularly  polarized  light.  The  strength  of  the  magnetic  field  determines  the 
magnitude  of  deflection  of  the  electrons  and  therefore  the  electron  z  velocity.  Varying 
the  electron  z  velocity  (i.e.  varying  the  electron  energy)  makes  an  FEL  tunable  over  a 
large  range  of  wavelengths.  The  undulator  parameter,  K  =  ©BXo/2nmc^,  where  E  is 
the  root-mean-square  (rms)  value  of  the  undulator  magnetic  field,  characterizes  the 
undulator’s  magnetic  field  strength.  In  a  linear  undulator  E  =  S/V2,  because  of  non¬ 
constant  acceleration,  while  in  a  helical  undulator  S  =  S,  because  of  the  continuous 
acceleration,  where  B  is  the  peak  magnitude  of  the  undulator  magnetic  field. 
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While  the  quantum  mechanical  approach  shows  that  the  emission  of  a  photon 
causes  a  recoil  in  the  electron,  a  classical  approach  to  the  FEL  theory  neglecting  recoil 
can  be  used  because  the  energy  loss  is  small.  The  ratio  of  energy  lost  with  electron 
initial  energy  is  on  the  order  of  10"^.  It  requires  many  such  emissions  to  create  an 
effect  that  would  be  significant,  due  to  the  high  energy  of  the  electrons.  With  the 
optical  beam,  the  photon  density  is  so  great  that  the  statistical  fluctuations  of  the 
optical  beam  can  also  be  ignored  [11]. 

The  classical  approach  describes  an  electron  wiggling  and  emitting  light  while 
transversing  the  undulator’s  magnetic  fields.  The  light  is  emitted  along  the  electron 
beam  axis  into  a  forward  cone.  Figure  2.1  [4]  shows  the  interaction  between  electrons 
and  the  optical  field.  It  shows  that  the  optical  electric  field  acting  on  an  electron  will 
take  energy  away  from  an  electron.  To  accomplish  this  transfer  of  energy  the  optical 
field  and  the  electron  must  have  the  correct  relative  phase.  The  spread  of  electron 
phases  in  comparison  to  the  optical  field  causes  some  of  the  electrons  to  take  energy 
from  the  optical  field  when  they  are  out  of  phase.  The  electrons  then  form  bunches 
over  the  optical  wavelength.  The  end  result  is  that  coherent  light  is  formed  from  these 
electron  bunches. 

As  mentioned,  unlike  classical  lasers  where  coherent  light  is  produced  by 
electrons  transitioning  from  one  quantum  state  to  another,  the  FEL  produces  light  by 
the  interaction  of  a  relativistic  electron  beam  and  external  forces.  The  FEL  is  then 
able  to  generate  coherent  light  at  tunable  frequencies  depending  on  the  undulator 
parameter  K.  As  the  electrons  are  relativistic  the,  undulator  wavelength  is  Lorentz 
contracted  to  X'o  =  Vy  in  fh©  electron  beam  frame.  The  optical  field  passing  over  the 
electrons  is  Doppler  shifted  to  a  wavelength  X'  =  (1  +  ~  2yX,  where  X  is  the 

optical  field  wavelength  in  the  lab  frame  [3].  To  have  sufficient  energy  exchange 
between  the  electron  and  the  optical  field  a  condition  of  resonance  must  occur, 
X'  =  Xfl'  in  the  beam  frame.  In  the  undulator  frame 

+  ^  (2.1) 
2f 

By  evaluating  Eq.  2.1,  it  is  apparent  that  an  FEL  can  be  continuously  tuned  by 
changing  either  the  electron  energy  ymc^,  the  undulator’s  magnetic  field  strength  B,  or 
the  undulators  wavelength  Xq-  The  undulator  can  also  be  designed  to  compensate  for 
the  falling  out  of  resonance  by  the  electrons  while  they  give  up  energy.  By  changing 
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Figure  2-1 .  Fundamental  physics  of  the  FEL  interaction. 


Iq  through  the  undulator  the  electron  can  be  kept  in  phase  with  the  optical  field  longer. 
This  is  called  tapering  the  undulator.  which  increases  the  efficiency  of  the  FEL. 

B.  ELECTRON  DYNAMICS  AND  THE  PENDULUM  EQUATION 

The  individual  electron’s  motion  inside  an  undulator  is  described  by  solving  the 
classical  equations  of  motion.  The  Lorentz  force  is  the  only  significant  force  acting  on 
an  electron  traveling  down  the  length  of  the  undulator  [3].  In  the  presence  of  the 
undulator  and  optical  fields,  the  Lorentz  force  equations  are 

^  =  ^[£',+px(g,  +  S„)l  (2-1a) 

at  me 

^  ^  .  e,)  (2.1b) 

at  me 

=  ’  (2.1c) 

where  ^  =  Vie  is  the  electron  normalized  velocity  vector.  The  magnetic  field  of  a 
helical  undulator  Sy  is  circularly  polarized  and  is  given  by 

%-B  [eosikoz),  sin(kQz),  0]  ,  (2.2) 

where  the  undulator  wave  number  ko  =  27i/Xo,  and  z  is  the  direction  of  the  undulator 
axis.  The  optical  electric  and  magnetic  fields  are  assumed  to  be  circularly-polarized 
plane  waves  and  are  given  by 

=  E  [cos(\y),  -s/n(\|/),  0]  ,  E  [s/n(\|;),  cos(\|/),  0]  .  (2.3) 

where  E  is  the  optical  field  magnitude  and  \j/  =  /cz  +  (of  +  <|).  The  optical  wave  number 
\s  k  =  QnIK,  the  optical  frequency  is  ©  =  ek,  and  ^  is  the  phase  of  the  optical  wave. 

Substituting  the  undulator  magnetic  field  Eq.  2.2  and  the  optical  wave 
electromagnetic  fields  Eq.  2.3  into  the  Lorentz  force  equations  Eq.  2.1  ,  they  become: 

4-('>^i)  =  —  [E(2-p^)(cos(M/).-s/n(v).0)  +  Pz  B{-sin{koZ),  eos(kQZ),  0)](2.4a) 
at  ^  me 

4(y^z)  =  —  [E(PxCOS(v)  -  pyS/n(M/))  -I-  B  {^^sin{koZ),  -  pyeos{koz)]  (2.4b) 
dt  me 

—  E  [PxCOs(m/)  -  PySin(M;)]  .  (2.4c) 

dt  me 

where  =  (Py,  py,  0).  As  the  electrons  travel  down  the  undulator  in  the  z  direction 
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they  are  acted  upon  by  a  transverse  force  from  the  optical  electric  and  magnetic  fields 
and  the  magnetic  field  from  the  undulator.  The  force  from  the  optical  field  is 
proportional  to  £  (1  -  Pz)  while  the  force  from  the  undulator  magnetic  field  is 
proportional  to  B.  For  a  relativistic  electron  p^  ^  1 ,  so  1  -  p^  — >  0.  Therefore,  we 
can  assume  E  (1  -  p^)  <  S  p^.  The  force  acting  on  the  electron  in  the  transverse 
direction  is  mainly  due  to  the  undulator  magnetic  field.  Equation  2.4a  can  then  be 
approximated  as 

^  Pz  B[-sin{koZ),  cos{kQz),  0]  .  (2.5) 

at  ^  me 


On  integrating  Eq.  2.5,  the  electron’s  velocity  in  the  transverse  direction  is  found  to  be 


[cos (/foZ),  sin(/foz),  0]  , 


(2.6) 


where  K  =  eB7^2nmc^  =  eBlkomc^  is  the  undulator  parameter. 

Examining  the  transverse  and  longitudinal  motion  of  the  electron  as  it  travels 
down  the  undulator,  it  is  noted  that  the  transverse  motion  is  small,  Pj^ «  K/y  1 . 

Inserting  Eq.  2.7  into  Eq.  2.5c  gives  the  electron  energy  evolution. 


“Y-  -  ■  Sr)  =  . 


dt 


me 


yme 


(2.7) 


where  the  dimensionless  electron  phase  in  the  combined  undulator  and  optical  fields  is 
t^  +  ^  =  \\r  +  kQ2  =  {k  +  ko)z  -  ©f  +  (|>.  The  initial  value  of  the  electron  phase  is 
^0)  =  Co  =  (^  +  kQ)zo.  Since  k  »/co  then  Co  =  2n^o/J^,  when  y »  1 .  The 
electron  phase  C  determines  the  electron’s  z  position  relative  to  a  single  optical 
wavelength  X.  By  using  Eq.  2.1c  and  inserting  in  Eq.  2.6  gives 

,  (2.8) 


or 

(1  +/<2)y-2=1  -P^2  . 

Then,  take  the  time  derivative  of  both  sides 

-2y-^Y(1+K2)  =  -2p^P^  , 

which  relates  y  to  p^  as 

Y  _  Y^PzPz 

Y  (l+K^)  ■ 


(2.9) 

(2.10) 

(2.11) 
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Taking  the  first  derivative  of  the  electron  phase  C  gives 


t  =  {k  +  ko)i  -  (0  =  (/f  +  /fo)cPz  -  ®  . 


(2.12) 


and  a  second  derivative  gives 


C  =  (k  +  /Co)cp, 


By  using  Eq.  2.11,  and  Eq.  2.13,  and  the  resonance  condition 

1  =  OC  _ _  oe  -L. 

Y  (1  +  K^{k  +  ko)  (o(1  +  K^)  2(00  ’ 


(2.13) 


(2.14) 


where  coq  =  co(1  +  K^)/2y^.  The  electron’s  equation  of  motion  Eq.  2.7  can  be  combined 
with  Eq.  2.14  to  give 

"  ^(JHqQKE  k  ^  w-\ 

'r  «  .  X\  /O  H  C\ 


C  =  ^ - COS(C  +  <|)) 

"fmc 


(2.15) 


Define  the  dimensionless  time  x  =  ctlL  where  L  is  the  length  of  the  undulator  to 
make  the  equation  of  motion,  Eq.  2.15,  dimensionless.  As  an  electron  passes  through 
the  undulator,  t  =  0-^  1.  The  dimensionless  optical  field  is  defined  as  a  =  lale'*, 
where  la  I  =  ^nNeKLEl-fmc^,  and  N  is  the  number  of  periods  in  the  undulator.  The 
pendulum  equation  can  be  written  as 

00  Q 

^  =v  =  la  I  cos(C +  <|))  .  (2.16) 

00  0 
where  ^  is  the  second  derivative  with  respect  to  x,  v  is  the  first  derivative  with  respect 

to  X,  and  the  dimensionless  electron  phase  velocity  is  defined  as 

V  =  L[{k  +  /co)p^  -  k].  When  v  =  0,  the  electrons  are  at  resonance  with  the  optical 

wave  and  A,  =  Xo(1  +  K^t^).  Changes  In  the  electrons’  energy  correspond  to  giving 

or  taking  energy  from  the  optical  wave.  As  the  electron  decreases  its  energy,  its 

phase  velocity  decreases  and  it  gives  up  energy  to  the  optical  wave.  The  maximum 

energy  loss  of  the  electron  occurs  when  ^  + 1|>  =  n. 


C.  THE  FREE  ELECTRON  LASER  WAVE  EQUATION 

The  development  of  a  complex  optical  wave  equation  starts  from  Maxwell’s  wave 
equation  acting  on  a  vector  potential  ^{r,t)  driven  by  the  current  density 

C  ot  ^ 
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where  +  3^/3/^  +  3^/32^,  the  vector  potential  ^(z,f)  is 


X(z,0  =  —  E{z,t)  [sin(/cz  -  cot  +  (|>(z.f)).  cos(/cz  -  cot  +  <l)(z.f)).  0]  .  (2.18) 

CO 


and  the  electric  field  ^  is 

^  =  E(z,f)[cos(/(z  -  cot  +  <t>(z,0),  -sin(/cz  -  cot  +  ^{z,t)),  0]  , 

where 

c  3f  ■ 


(2.19) 


(2.20) 


Since  the  wave  is  traveling  only  in  the  z  direction  the  Maxwell  wave  equation 
operator  then  is 

- ^  .  (2.21) 

3z^  3f^ 

Taking  the  first  derivative  of  the  vector  potential  X  with  respect  to  z  gives 

[cos(m/),  -sin(\y),  0]  .  (2.22) 


M.lM[si„(^),CO8(v),0]+| 


k+^ 

Bz 


where  \\f  =  kz  -  (at  +  ^,  and  c/©  =  Mk.  Taking  the  second  derivative  with  respect  to  z 
gives 

3^^  1  BE 


Bz^  k  Bz 


k  +  ^ 

dz 


* 


[cos(v).  -sin(M/).  0] 

[cos(m/).  -sin(\(f),  0] 


k  +  ^ 

dz 


+  [cos(\|/),  -sin(i|/),  0]  +  -f  (/f  +  [-sin(\j/),  -cos(\|/),  0]  .  (2.23) 

Assuming  slowly  varying  amplitudes  and  phases  where  BElBz  «  kE,  B^lBz  <«:  k^, 
BElBt  <  (oE,  B^lBt  <  co<t>.  and  ©  =  /fc  allows  us  to  use  only  first-order  terms  and 
simplify  Eq.  (2.23)  to 


3^^  ^  o  3£ 
Bz^  ~  dz 


2  «  2  ^[cos{\v).  -sin(M;),  0]  + 


k  +  ^ 

Bz 


[-sin(\|/),  -cos(\j/),  0]  .  (2.24) 


By  expanding 


k  +  ^ 

k^  +  2  k^  + 

3(|) 

Bz 

Bz 

dz 

(2.25) 
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knowing  (3<|>/0z)^  <  d^tdz.  and  putting  Eq.  2.25  into  Eq.  2.24,  it  becomes 

[-sin(\i/),  -cos(\)/),  0]  .(2.26) 


[cos(m/),  -sin(M/),  0]  + 

dz^  oz 


k  E  +  2E  ^ 
dz 


The  time  derivative  portion  of  the  Maxwell’s  wave  equation  operator  is 

-E  [sin(\|r),  cos(\4;),  0]  .  (2.27) 


1  c 


dt^ 


CO 


Taking  the  first  derivative  with  respect  to  time  gives 

[E  (sin(M/),  cos{\\f),  0)]  =  ■^[sin(\|/),  cos(Mr),  0] 


+  E 

The  second  derivative  gives 

-|^[sin(vi/),  cos(\|/),  0]  = 


a<i)  _ 
dt 


CO 


[cos(vi/),  -sin(\^/),  0]  . 


(2.28) 


a<i>  _ 
dt 


CO 


^2p 

cos(mc),  -sin(M/),  0]  +  — ^[sin(\|/),  cos(vj/),  0] 
at 


dt 


a<i)  _ 
dt 


(O 


[cos(Mr),  -sin(\K),  0]  +  E^[cos(\j/),  -sin(\i;),  0] 

dt 


+  E 


a(i>  _ 
dt 


CO 


[-sin(\|/),  -cos(m/),  0] 


(2.29) 


Again,  assuming  the  slowly  varying  amplitudes  and  phases  allows  us  to  neglect  all  but 
first-order  terms,  so  Eq.  2.29  becomes 

12 


«  -2  c  ■^[cos(\y).  -sin(\|;),  0]  +  -^  E 
dr  ar  CO 


a<i>  _ 

dt 


CO 


[-sin(M;),  -cos(\|/),  0]  (2.30) 


By  expancfing 


-^-co 

2 

a<i» 

dt 

dt 

w  J 

+  ^-^2 
dt 


(2.31) 


again  knowing  (d^/dt)^  <  d^ldt,  we  put  Eq.  2.31  into  Eq.  2.30  and  it  becomes 

^  js  -2  c  ^[cos{\v).  -sin(\|/),  0]  +  c  (co  E  -  2  E  ^)[-sin(\j/),  -cos(\k),  0]  (2.32) 


de  at 

Putting  Eq.  2.26  and  Eq.  2.32  into  Eq.  2.17  gives 
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tz  ,  (2.33a) 


d^  1 

CM 

li 

1^  +  1 

dE 

ti  +  2  E 

+  1 

dz^  dt^ 

dz  C 

dt 

dz  C  dt 

where  the  unit  vectors  are  =  [cos(\}/),  sin(M/),  0]  and  ^  =  [-sin(\(/),  -cos(m/),  0].  We 
also  know  Eq.  2.33a  is  a  multiple  of  the  current  density  in  Eq.  2.17  so 


471 


2 


^  1  3E 

dz  C  dt 


ti  +  2  E 


dz  C  dt 


(2.33b) 


To  solve  for  the  time  variation  of  the  field  envelope  make  the  change  of  variables 
z  ^  2  -  c  t.  Then  2  =  z  +  c  t  where  the  undulator  length  L  =  NXq,  N  is  the  number 
of  periods,  Ao  is  the  wavelength  of  the  undulator,  and  t  =  ct/L  is  the  dimensionless 
time.  By  making  the  change  of  variables,  we  change  the  partial  derivatives, 

A  =  .^A  +  1lA=iA  +  o  ,  (2.34) 

dz  dz  02  dz  dx  02 


and 


C  dt  C  dt  dz  C  dt  dx  02  L  dx  ' 
Putting  Eq.  2.34  and  Eq.  2.35  together  simplifies  the  differential  operator  to 

_L  +  1  =  l  A 

dz  C  dt  L  dx  ' 

Making  the  substitution  Eq.  2.36  into  Eq.  2.33b  and  projecting  onto  gives 
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i  M 

L  dx  ■ 


(2.35) 


(2.36) 


(2.37) 


Doing  the  same  to  Eq.  2.33b  but  projecting  onto  tz  gives 

-^J^_-tz^2E  .  (2.38) 

The  current  density  for  a  single  electron  is  =  -  The  electron  trajectory 
will  be  found  to  be  ^4^  =  -  {Kly)  [cos(/foz),  sin(kor),  0],  so  that 

;t  *  [cos(/foz),  sin(/foz).  0]  ,  (2.39) 

+  y 

where  y=  I/V1  -  v^/c^,  v  is  the  velocity  of  the  electron,  /cq  is  the  undulator  wave 
number  2ti/Xq,  and  the  undulator  parameter  is  K  =  eE/mkoC  . 
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Inserting  Eq.  2.39  and  the  value  for^i  into  Eq.  2.37  gives 

M  _  2ti/^  [cos(/foz).  sin(/coz),  0]  •  [cos(m/),  sin(\y).  0]  .  (2.40) 

dx  Y 

On  performing  the  dot  product,  Eq.  2.40  becomes 

M  =  _  +  (2.41) 

dx  y 

where  the  electron  phase  is  C  =  {k  +  l<o)z  -cof- 

The  wave  equation  Eq.  2.41  is  for  a  single  electron.  The  average  value  over 
many  electrons  gives 

^  =  -  <cos(C  +  (|>)>  ,  (2.42) 

dx  Y 


where  p  is  the  electron  particle  density  and 

Ne 

<cos(C  +  <!>)>=£ 

/=0 


cos(C,  +  ^i) 
Ne 


(2.43) 


Similar  manipulations  of  Eq.  2.38  gives  us 

^  <sin(C  +  ({.)>  .  (2.44) 

dx  Y 

The  complex  electric  field  is  E  =  lEl  e'  and 

M  =  Ml.  e'^  +  i  \E\  ^  .  (2.45) 

8x  dx  dx 

When  Eq.  2.42  and  Eq.  2.44  are  inserted  into  Eq.  2.45,  we  get 

M  =  _  2Ti/.e/Cp  (<gQ5(^  +  _  /  <sin(C  +  <j))>]  e'  ,  (2.46) 

dx  Y 

which  is  also 

d^  __  2nLeKp  ,{>  (2.47) 

dx  Y 

From  the  pendulum  equation,  we  know  that  our  dimensionless  optical  wave  field 
is  a  =  la  I  e'  ^  where  la  I  =  4TiA/KLe  I E  l/y^ mc^,  so  that 

^  a  =  Sv^e^pNK^L^  (C4<|.)>  <|.  (2.48) 

9x  Y^mc^ 
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The  dimensionless  current  is  j  =  Sv^e^NK^L^I'^mc^.  Which  leads  us  to  the 
Free  Electron  Laser  optical  wave  equation 

a  =  -  y  <e-'  .  (2.49) 

Looking  at  the  real  and  imaginary  parts  of  Eq.  2.49,  you  have 

la  I  =  -j  <cos(C  +  (|))>,  I  =  -r^  <sin(C  +  <|>)>  ,  (2.50) 

lai 

where  la  I  is  the  optical  field  amplitude  and  ^  is  the  optical  field  phase.  Equation  Eq. 
2.50  shows  that  bunching  of  electrons  in  phase  at  (C  +  <|>)  =  n  drives  the  optical  wave 
amplitude  and  produces  a  gain,  while  bunching  electrons  around  (C  +  (!>)  =  n/2  will 
drive  the  optical  phase.  Increasing  the  electron  beam  current  density  j  will  increase 
gain  and  optical  phase  evolution.  The  optical  fields  increase  in  magnitude  will  reduce 
the  effects  of  the  electron  phase  evolution. 


D.  LOW  GAIN  FEL 

The  Maxwell-Lorentz  (standard  classical,  i.e.  non-quantum  dynamics)  theory  for 
an  FEL  can  be  formed  by  combining  the  simple  pendulum  equation,  Eq.  2.16,  and  the 
optical  wave  equation,  Eq.  2.49.  The  strength  of  the  dimensionless  current  density  j 
governs  the  coupling  between  electrons  and  light.  It  determines  the  response  of  the 
optical  wave  to  the  bunching  of  electrons  in  phase  C-  The  Maxwell-Lorentz  theory  is 
generally  valid  for  both  weak  (la  I  <«:  n)  or  strong  (la  I  >  n)  optical  fields  [5]. 

From  the  pendulum  equation,  Eq.  2.16,  it  is  evident  that  the  electron  energy  loss 
is  greatest  when  the  electron  bunch  is  at  ^  +  <|)  =  n,  which  provides  for  the  maximum 
gain  in  the  optical  field.  As  the  electron  bunch  moves  toward  C  +  <l>  =  0,  over-bunching 
occurs  and  energy  is  absorbed  back  into  the  electron  beam. 

Maximum  coupling  between  the  optical  field  and  the  electron  beam  occurs  when 
the  initial  electron  phase  velocity  is  at  resonance,  vq  =  0  [6].  A  "realistic"  beam  of 
monoenergetic  electrons  will  enter  the  undulator  with  a  random  spread  in  phase  C- 
There  will  be  just  as  many  electrons  gaining  energy  from,  as  losing  energy  to,  the 
optical  field,  resulting  in  zero  gain.  To  achieve  a  useful  gain,  the  electrons  must  have 
a  velocity  slightly  above  resonance  so  that  more  electrons  give  up  energy  than  absorb 
it  in  an  exchange  with  the  optical  field. 
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To  understand  the  operatic  ?  of  an  FEL,  it  is  useful  to  examine  the  evolution  of 
the  electrons  as  they  move  through  phase  space  (C,v).  The  path  in  phase  space  for 
each  electron  is  dependent  on  the  initial  values  of  Co  ^nd  vg.  A  low-gain  FEL 
produces  an  electron  phase-space  path  that  can  be  represented  by  a  simple 
mechanical  pendulum  [3],  and  is  shown  analytically  by 

=  vg^  -  2  la  I  [sin(C  +  <|>)  -  sin(Cg)]  .  (2.51) 

By  examining  Eq.  2.51 ,  we  can  note  that  there  are  fixed  points,  (n/2,  0)  and  {37c/2,  0), 
where  the  pendulum  does  not  evolve.  The  electrons  near  these  phase  space  points 
evolve  slowly.  The  electrons  that  are  at  the  stable  fixed  point  (7x/2,  0)  correspond  to  a 
pendulum  at  the  bottom  of  its  arc.  The  electrons  with  a  large  phase  velocity 
Ivl  >  2lal^'^  are  in  an  open  orbit.  This  would  correspond  to  a  pendulum  swinging  in 
one  direction  continuously  about  its  pivot  point.  An  electron  in  a  closed  orbit  is  like  a 
pendulum  swinging  with  a  small  amplitude  oscillating  about  the  stable  fixed  point.  A 
change  in  (|>  would  cause  a  change  in  the  horizontal  position  of  the  pendulums  fulcrum. 
The  phase-space  path  that  separates  the  open  and  closed  orbits  of  an  electron  is 
called  the  "separatrix"  and  is  given  by 

=  21a  l[1-sin(Cs +<!>)]  -  (2.52) 

The  separatrix  has  a  peak-to-peak  height  of  4la  1^^^.  Its  horizontal  position  depends  on 
the  optical  phase  (|i.  As  the  optical  field  grows,  so  does  the  separatrix  height. 

Since  the  electron’s  evolution  in  phase  space  is  periodic  over  each  optical 
wavelength  X,  the  interaction  in  an  FEL  can  be  modeled  by  tracking  sample  electrons 
over  one  wavelength.  Figure  2-2  is  an  example  of  a  phase  space  plot  with  low  current 
y  =  1,  a  moderate  optical  field  strength  of  lal  =  4,  and  an  electron  phase  velocity  of 
Vo  =  2.6. 

The  electrons  enter  the  undulator  at  x  =  0,  in  a  uniform  distribution  of  initial 
phases,  Co  =  ~  ^2  3n/2.  The  electrons  then  interact  with  the  optical  field  of 
strength  ag  =  4.  The  initial  electrons  are  displayed  as  light  gray  dots.  As  they  evolve 
through  the  undulator  they  become  darker  until  they  become  black  dots  at  the  end  of 
the  undulator,  t  =  1 .  The  electrons  to  the  left  experience  an  increase  in  phase  velocity 
which  causes  them  to  move  to  the  right,  while  the  electrons  to  the  right  experience  a 
loss  of  energy  and  decrease  in  phase  velocity.  This  produces  a  spatial  bunching 
observed  near  phase  =  n  at  x  =  1 ,  so  more  electrons  are  giving  up  energy  than 
taking  energy  away  from  the  optical  field.  There  is  a  limit  in  the  process  of  energy 
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Figure  2-2.  Phase-space  plot  for  electrons  entering  with  vq  =  2.6. 


exchange  where  the  electrons  have  given  up  so  much  energy  that  their  phase  and 
phase  velocity  cause  the  electrons  to  start  absorbing  energy  back  from  the  optical 
field.  An  undulator’s  length  can  be  designed  to  prevent  the  electrons  from  reaching 
this  stage. 

The  optical  gain  is  defined  as  the  fractional  increase  in  optical  field  power  over 
time  X, 


G(t)  = 


(a(x)2  -  ao^) 


(2.53) 


Gain  is  pictorially  cfisplayed  in  the  upper  right  graph  of  the  phase  space  Figure  2-2. 
The  electrons  in  the  beam  are  assumed  to  be  uniformly  cfistributed  in  phase  Co 
have  the  same  initial  phase  velocity  vq.  A  change  in  electron  energy  corresponds  to  a 
change  in  electron  phase  velocity.  The  average  energy  lost  by  an  electron  is 
ymc^<v>  -  vo)/4nA/.  The  energy  lost  by  the  electron  determines  the  optical  gain  to  be 
G  =  2y(vo  -  <v>)la^  [3].  The  phase  velocity  average  is  found  by  expanding  the 
pendulum  equation,  Eq.  2.16.  The  small  signal  gain  equation  becomes  [3] 
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If  the  phase  velocity  is  far  aw  y  from  resonance,  IvqI  the  final  gain  at  x  =  1  will  be 
small.  Optimum  coupling,  resulting  in  maximum  gain,  occurs  when  IvqI  S  Ji-  The 
natural  gain  bandwidth  for  an  FEL  is  derived  from  the  relation  Avq  «  4-nNAyly.  The 
natural  gain  bandwidth  is 


At  ~J_ 

y  ~  2N  ' 


(2.55) 


An  FEL  gain  spectrum  and  optical  phase  shift  is  plotted  in  Figure  2-3,  for  a  weak 
field  ao  =  i,  low  current  y  =  1 ,  and  the  initial  phase  velocity  between  the  interval 
-1 2  ^vq  ^  1 2.  The  gain  curve  is  antisymmetric  about  vq  with  maximum  gain  at 

vo  =  2.6. 


Gain  and  Phase  Curves 

j=l  a  =1  N=25 

o 

0.135 

0.0 

-0.135 
0.08 

0,0 

-0.08 

Figure  2-3.  Weak-field,  low-gain  FEL  gain  spectrum  and  optical  phase  shift  plot. 
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E.  ENERGY  SPREAD 


The  previous  section  dealt  with  a  perfect  beam  where  the  electrons  were 
uniformly  distributed  over  the  optical  wavelength  with  the  same  initial  energy.  This  is 
not  achievable,  since  there  is  a  trade-off  between  electron  beam  quality  and  current 
density.  Increasing  current  /  tends  to  decrease  beam  quality  from  the  accelerator, 
while  increasing  the  length  of  the  undulator  by  increasing  N  increases  the  FEL’s 
sensitivity  to  beam  quality  because  the  natural  gain  bandwidth  is  narrower.  Where  the 
natural  gain  bandwidth  is  defined  as  the  range  of  phase  velocities  where  the  growth 
rate  is  reduced  by  10  %  below  its  peak  value  [3]. 

There  are  two  qualities  that  describe  electron  beam  quality:  the  electron  beam 
energy  spread,  Ay/y,  and  emittance,  e.  Electrons  with  same  initial  phase  C.  but 
different  phase  velocities  will  drift  apart  on  the  z  axis  [6].  At  the  end  of  the  undulator 
their  final  phase  difference  is  roughly  AC  «  Avt,  where  t  =  1 ,  for  each  electron.  It  is 
apparent,  for  a  change  in  initial  energy  Aymc^,  there  is  a  corresponding  change  in 
electron  phase  velocity,  Av  «  AnNAyly.  A  Gaussian  distribution  of  energies  with  an 
rms  spread  ,  Ay/y,  is  characterized  by  standard  deviation  oq  =  AnNAyly.  When 
gq  ~  Jt,  bunching  becomes  impaired  [6].  Typical  energy  spreads  from  an  accelerator 
are  on  the  order  of  Ay/y  «  0.001 . 

Emittance,  e,  is  the  average  angular  spread  times  the  average  beam  radius, 
e  =  TQ,  where  r  is  the  average  beam  radius  and  0  is  the  average  angular  divergence 
of  the  electron  beam.  The  normalized  emittance,  £„  =  ey  =  r0y,  is  a  measure  for  beam 
quality.  As  the  electrons  are  accelerated  and  beam  energy  increases,  the  normalized 
emittance  remains  constant.  There  is  a  Gaussian  spread  in  angles  with  rms  value  0 
resulting  in  an  exponential  distribution  of  phase  velocities.  The  standard  deviation  in  0 
is  Oq  =  AnN'f^lii  +  K^).  As  with  gq,  when  oq  is  equal  to  or  greater  than  n.  electron 
bunching  is  impaired. 


F.  HIGH-CURRENT  HIGH-GAIN  FEL 

In  high  current,  high  gain  FELs  j  »  n,  the  magnitude  of  the  optical  field  changes 
significantly  over  a  single  pass  of  the  electrons  through  the  undulator.  This  is  evident 
in  the  wave  equation,  Eq.  2.48.  As  we  did  in  the  low  gain  FEL.  assume  a  perfect 
quality  beam  operating  on  resonance,  vo  =  0.  From  the  pendulum  and  wave 
equations,  the  optical  field  grows  exponentially  for  large  j  as  it  travels  down  the 
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undulator  [3]: 


la(^)l  «  fo  e(y72)^N3x/2  . 


and  the  phase  evolves  as: 


The  FEL  gain  for  large  j  is 

G(^)  «  le072)’'^^2 

3 


(2.56) 


(2.57) 


(2.58) 


For  an  electron  beam  that  is  off  resonance,  the  single  pass  gain  is 

=  .  (2.59) 

9 

where  xg  =  (2/;)^^  is  the  bunching  time.  The  peak  gain  will  occur  slightly  above 
resonance.  The  phase  velocity  has  now  become  a  function  of  current.  As  y  ->  «>  the 
gain  spectrum  becomes  symmetric  centered  at  resonance. 

The  gain  bandwidth  is  Avg  *  The  bandwidth  for  a  low  gain  FEL  is 
Avfl  «  nl2.  It  is  evident  that  the  bandwidth  for  the  high  gain  case  can  be  much  larger 
than  that  for  the  low  gain  FEL. 


G.  THE  KLYSTRON  UNDULATOR 

The  purpose  of  a  "dispersive  section"  in  a  klystron  undulator  is  to  increase 
optical  gain  in  a  weak  optical  field.  There  are  three  major  sections  to  a  klystron 
undulator.  The  undulator  that  was  described  earlier  is  divided  into  two  sections 
separated  by  a  "drift  section"  or  a  "dispersive  section".  The  first  section  of  the 
undulator  is  called  the  "modulator".  It  develops  the  phase  velocity  differences  to 
prepare  the  electrons  for  bunching.  The  next  portion  is  either  a  drift  section,  i.e.  an 
open  section  that  allows  the  electrons  to  continue  their  bunching  without  interaction 
with  an  optical  field,  or  a  dispersive  section,  made  up  of  magnets  that  force  the 
electrons  far  off  resonance.  The  final  section  of  the  undulator  is  called  the  "radiator", 
where  the  electrons  and  the  optical  field  interact  and  create  the  coherent  optical  beam. 
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Mathematically,  there  is  no  difference  between  the  use  of  a  drift  section  or  a 
dispersive  section  [3].  The  drift  section  requires  a  large  distance  over  which  the 
electrons  travel  with  no  optical  interaction.  This  is  easily  achieved,  but  can  possibly 
double  the  length  of  the  undulator,  or  even  larger.  A  dispersive  section  forces  the 
electrons  off  centerline  in  an  arc,  and  allows  for  the  use  of  a  shorter  undulator.  The 
strength  of  either  type  of  system  is  given  by  the  dimensionless  drift  time  D.  where  D 
can  be  determined  by  the  number  of  undulator  periods  equivalent  to  the  drift  space. 

0  =  -^  .  (2.60) 

where  is  the  length  of  the  drift  space  in  undulator  periods. 

While  the  electrons  are  traveling  in  the  modulator  and  radiator  the  equations 
governing  their  motion  are  still  Eq.  2.16  and  2.49.  While  in  the  dispersive  or  drift 
section  there  is  no  electron-optical  wave  interaction.  The  only  change  in  electron 
phase  is  caused  by  the  constant  phase  velocity,  v;  therefore,  the  change  in  phase 
velocity  is  Av  =  0,  so  =  vD  is  the  only  change  in  electron  phase.  The  electrons 
then  will  bunch  without  any  interaction  from  the  optical  field,  and  the  optical  field 
strength  will  not  grow. 

When  using  a  klystron  the  gain  and  phase  equation  are  [3] 

C3(vo)  ®  sin(voD)  ,  (2.61) 

and 

<l»(vo)  «  ^  cos(voD)  .  (2.62) 

It  is  apparent  that  the  maximum  gain  is  about  jDl4  when  vq  =  rtl2D .  It  is  also 
apparent  that  with  low  j,  the  addition  of  a  klystron  enhances  gain.  The  spread  in  initial 
phase  velocity  is  critical.  The  natural  gain  bandwidth  becomes  narrower,  dictating  that 
beam  quality  improve.  The  spread  in  initial  phase  velocity  should  be  kept  to 
Avo  ^  which  is  smaller  than  that  required  for  a  regular  undulator.  The  stronger 
the  dispersive  section,  the  greater  the  need  for  increased  beam  quality. 

Figures  2-4a  and  2-4b  are  the  optical  gain  spectra  for  a  low-current  (y  =  1), 
weak-field  {Sq  =  1),  FEL  without  a  klystron  and  with  a  klystron,  D  =2.  Along  with  the 
narrowing  of  the  natural  bandwidth,  gain  is  significantly  improved  by  addition  of  a  drift 
section  or  dispersive  section. 
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Figures  2-5a  and  2-5b  are  the  optical  gain  spectra  for  a  high-current  (y  =  1 00), 
weak-field  (ao  =  1)  FEL  without  a  klystron  and  with  a  klystron  of  magnitude  D  =2. 
The  natural  bandwidth  is  increased  and  gain  is  significantly  improved  with  the  use  of  a 
klystron  undulator  even  in  a  high-current  FEL.  The  size  of  D  for  a  high-current  FEL  is 
significantly  smaller  than  that  for  a  low-current  FEL. 


****  Gain  Curve  **** 


-12  12 


Figure  2-4a.  Low-current,  weak-field  FEL  gain  spectrum,  without  dispersive  section. 


-0.572 


-12  12 
Figure  2-4b.  Low-current,  weak-field  FEL  gain  spectrum,  with  dispersive  section. 
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****  Gain  Curve  **** 

j=100  a  =1 
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75.67 

0.0 

-75.67 

Figure  2-5a.  High-current,  weak-field  FEL  gain  spectrum,  without  dispersive  section. 


****  Gain  Curve  **** 
j=100  aQ=l  D=2 


765.4 
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Figure  2-5b.  High-current,  weak-field  FEL  gain  spectrum,  with  dispersive  section. 


III.  SLAC  X  RAY  FREE  ELECTRON  LASER 


A.  INTRODUCTION 

Research  in  the  area  of  high  power  FEL’s  is  important  to  the  military  for  the 
development  of  missile  defense.  The  research  and  development  that  is  being 
performed  at  SLAC  is  in  the  high  peak  power  (GeV)  FEL  regime.  By  developing  hard 
X  ray  laser  light  at  a  wavelength  of  1.5  A,  imaging  of  DNA  base-pairs  may  be 
accomplished  [7].  The  FEL  is  the  best  possible  source  for  an  X  ray  laser. 

The  SLAC  project  proposes  to  develop  coherent  hard  X  ray  radiation  at 
wavelengths  as  short  as  1 .5  A.  It  will  use  a  single  pass  of  an  electron  beam  through 
an  undulator  with  no  initial  optical  field.  The  optical  field  will  be  started  from  noise 
using  a  process  called  Self-Amplified  Spontaneous  Emission  (SASE).  There  are  no 
mirrors  for  feedback  of  the  optical  beam.  The  electron  beam  makes  a  single  pass 
through  the  undulator  creating  a  cone  of  light  in  the  direction  of  the  electron  beam’s 
path.  The  only  interaction  between  the  electron  beam  and  the  optical  field  happens 
during  this  single  pass.  Figure  3-1  is  a  diagram  of  the  system  proposed  by  SLAC.  It 
combines  the  use  of  the  SLAC  linear  accelerator  (linac)  with  an  undulator  to  produce 
the  coherent  X  ray  light.  The  mirror  merely  directs  the  light  beam  at  grazing  angles. 


Rgure  3-1.  Schematic  of  the  SLAC  linac  driving  a  coherent  X  ray  FEL. 


The  system  design  would  incorporate  a  section  of  the  SLAC  linac  accelerating 
the  electron  beam  to  an  energy  of  ymc^  =  15  GeV  with  peak  current  of  /  =  5000  A. 
Table  3-1  contains  the  design  parameters  for  two  systems  that  would  produce  either  a 
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4.5  A  or  1 .5  A  laser  beam  [8]. 


PARAMETERS 

A 

B 

X  (Radiation  Wavelength) 

1.5  A 

4.5  A 

Xq  (Undulator  period) 

2.67  cm 

4.0  cm 

B  (Peak  undulator  field) 

1.173  T 

1.66  T 

Undulator  type 

Helical 

Linear 

K  Undulator  parameter  (peak) 

3.4 

4.2 

ymc^  (Electron  energy) 

15GeV 

15GeV 

/  (Peak  current) 

5000  A 

5000  A 

Ay/y  (Energy  spread) 

0.0002 

0.0002 

L  (Undulator  length) 

25  m 

38  m 

Output  power 

11  GW 

66  GW 

Table  3-1 .  Proposed  parameters  for  an  X  ray  PEL  at  SLAG. 


The  dimensionless  current  for  /  =  5000  A  is  j  =  2056  for  the  4.5  A  laser  and 
y=3121  for  the  1.5  A  laser.  As  mentioned  earlier,  the  dimensionless  current 
determines  the  response  of  the  optical  wave  to  the  electron  beam  bunching.  The  gain 
of  the  system,  because  of  the  high  y,  is  in  the  high-gain  regime.  Rgure  3-2  is  a 
phase-space  plot  for  the  4.5  A  system  and  figure  3-3  is  for  the  1 .5  A  system.  Each 
simulation  is  started  from  an  optical  field  ao  =  0  and  the  electrons  are  all  started  with 
random  initial  phase  C  between  -7t/2  and  3nl2.  The  simulations  ran  with  an  undulator 
length  of  L  =  35  m.  A  filling  factor  F  =  0.1  was  used,  where  the  filling  factor  is  the 
ratio  of  the  area  of  the  electron  beam  to  the  optical  beam.  The  electron  beam  radius 
Tf,  =0.0015  cm.  The  initial  electron  phase  velocity  has  a  standard  deviation  of 
oq  =  2.2  and  4.2,  where  oq  -  4nNAyly.  The  average  initial  phase  velocity  is  centered 
around  vq  =  0.  The  electrons,  in  both  scenarios,  have  bunched  and  are  on  the 
bottom-half  in  the  separatrix  at  the  end  of  the  undulator,  t  =  1 .  The  electrons  have 
given  energy  to  the  optical  field.  The  phase-space  simulations  indicate  the  production 
of  an  X-ray  beam  with  saturation  of  the  system  at  undulator  lengths  of  34  and  35 
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meters,  respectively.  The  logarithmic  power  evolution  shown  in  the  upper  right-hand 
corners  of  each  figure  indicate  the  systems  have  reached  saturation.  The  peak  power 
achieved  for  the  4.5  A  is  approximately  66  GW  and  for  the  1 .5  A  it  is  approximately  1 1 
GW. 


Figure  3-2.  Phase-space  plot  for  4.5  A  system. 


Figure  3-3.  Phase-space  plot  for  1.5  A  system. 
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In  either  seer  1o,  the  proposed  parameters  would  create  a  hard  X  ray  FEL  as 
indicated  by  the  pi;  se-space  simulations.  A  large  portion  of  the  FEL  cost  is  in  the 
production  of  the  undulator;  reducing  the  length  of  the  undulator  reduces  the  FEL  cost. 
A  method  of  reducing  the  length  of  the  undulator  is  by  introducing  a  dispersive  section 
in  the  undulator  to  make  an  FEL  klystron.  A  penalty  resulting  from  the  klystron  design 
is  the  onset  of  saturation  at  weaker  fields  than  without  the  klystron  [9]. 


B.  KLYSTRONS  AND  HIGH-CURRENT  HIGH-POWER  FEL 


An  FEL  with  an  optical  wavelength  of  A,  =  40  A  was  proposed  uang  the  SLAG 
linac  as  an  electron  beam  source  [10].  Table  3-2  gives  tie  original  proposed 
parameters  for  the  optical  wavelength  of  X  =  40  A.  The  system  required  an  undulator 
with  a  length  of  60  meters.  The  phase-space  plot  in  Fig.  3-4  shows  a  simulation  of  a 
soft  X-ray  FEL  at  40  A  using  the  parameters  in  Table  3-2,  an  initial  radiation  field 
ao  =  0,  and  a  phase  velocity  spread  of  standard  deviation  oq  =  3.6  centered  around 

Vo  =  0. 


PARAMETERS 

VALUE 

X  (Radiation  Wavelength) 

40  A 

Xo  (Undulator  period ) 

8.3  cm 

K  (Undulator  parameter) 

4.2 

Undulator  type 

Linear 

En  (Normalized  rms  emittance) 

3.0  mm-mrad 

yme^  (Electron  energy) 

7.0  GeV 

/  (Peak  current) 

2500  A 

Energy  spread 

0.0004 

L  (Undulator  length) 

60  m 

Output  power 

28  GW 

Table  3-2.  Proposed  parameters  for  a  X  =  40  A  X  ray  FEL  at  SLAG. 
(Without  dispersive  section) 
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***  FEL  Phase  Space  Evolution  *** 
j=4000  ^0=0  ''o~®  H=732 


Figure  3-4.  Ph£ise-space  plot  for  40  A  system. 
(Without  dispersive  section) 


The  optical  field  grows  from  noise  as  a  result  of  interaction  with  200,000 
electrons  randomly  positioned  in  the  initial  electron  phase.  The  dimensionless  current 
j  =  4000  corresponds  to  a  current  /  =  2500  A.  Notice,  at  the  end  of  the  undulator, 
X  =  1 ,  the  electrons  have  already  reached  saturation  as  shown  in  the  power  plot  in  the 
upper  right  diagram.  Having  gone  through  one  complete  transversal  inside  the 
separatrix,  the  bunched  electrons  are  near  the  bottom  of  the  separatrix,  for  a  second 
time,  where  maximum  energy  transfer  to  the  optical  beam  occurs.  When  the  electrons 
are  bunched  at  the  upper  half  of  the  separatrix  the  electron  beam  has  taken  back 
energy  from  the  optical  beam  as  indicated  by  the  drop  of  the  optical  power  after  the 
onset  of  saturation.  Reducing  the  undulator  length  to  approximately  40  m,  at  the  onset 
of  saturation,  will  result  in  peak  power  at  the  end  of  the  undulator.  The  majority  of 
electrons  will  be  at  the  bottom  of  the  separatrix  and  the  electrons  will  have  given  up 
the  maximum  energy  possible  to  the  optical  beam.  By  reducing  the  length  of  the 
undulator  we  reduce  the  dimensionless  current,  j  =  Sn^e^pNK^L^Fl'fmc^.  With  an 
undulator  length  of  40  m,  the  dimensionless  current  is  j  -  2000.  Figure  3-5  is  the 
phase-space  plot  for  such  a  system. 

To  reduce  costs  of  manufacturing  an  undulator  of  such  a  long  length  and 
maintaining  the  same  results,  a  dispersive  section  can  be  added  to  an  undulator,  thus 
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making  a  klystron.  Table  3-3  lists  the  proposals  for  a  system  that  uses  a  klystron  to 
produce  a  40  A  PEL  [1 1], 


***  FEL  Phase  Space  Evolution  *** 

j=2000  a  =0  v„=0  N=723 

o  o 


00=3 . 6 


-n/2  C  37t/2  0x1 


Figure  3-5.  Phase-space  plot  for  40  A  system,  undulator  length  of  40  m. 

The  FEL  design  at  A.  =  40  A  with  a  klystron  has  an  undulator  length  of  L  =  17 
meters,  electron  beam  energy  of  5  GeV,  and  a  peak  current  of  5000  A,  which 
corresponds  to  a  dimensionless  current  density  of  j  =  932.  This  FEL  was  examined 
without  a  klystron  and  then  with  a  klystron. 

The  case  without  the  klystron  is  shown  in  Fig.  3-6.  The  electron  bunch  is  at  a 
phase  where  it  gives  energy  to  the  optical  field  at  the  end  of  the  undulator.  The 
dimensionless  power  P  corresponds  to  a  peak  power  of  1 5  GW.  In  Rg.  3-7,  the  use 
of  the  klystron  FEL  with  the  same  parameters  and  D  =  0.5  shows  that  electrons  have 
bunched  in  the  lower  half  of  the  separatrix  and  at  the  optimal  phase  relation  ^  =  n 

for  maximum  gain.  The  peak  power  calculated  for  the  klystron  FEL  is  31  GW.  While 
not  shown,  saturation  will  occur  earlier  with  the  use  of  a  dispersive  section  than 
without,  reducing  the  overall  length  of  the  undulator. 
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PARAMETERS 

VALUE 

A.  (Fund.  Rad.  Wavelength) 

40  A 

Aq  (Undulator  period) 

4.0  cm 

K  (Undulator  parameter) 

4.2 

Undulator  type 

Linear 

c„  (Normalized  rms  emittance) 

3.0  mm-mrad 

ymc^  (Electron  energy) 

5.0  GeV 

/  (Peak  current) 

5000  A 

Energy  spread 

0.0004 

L  Undulator  length 

17m 

Output  power 

2  GW 

Table  3-3.  Proposed  parameters  for  a  A,  =  40  A  X  ray  PEL, 
(With  dispersive  section). 


***  FEL  Phase  Space  Evolution  *** 
j=932  v^=0  N=413 
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Rgure  3-6.  Phase-space  plot  for  40  A  X  ray  FEL, 
(Without  dispersive  section). 
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To  generalize  the  optical  field  in  the  longitudinal  direction,  several  sites  along  the 
complex  wave  envelope,  a(z),  are  followed.  The  length  in  the  z  direction  is  scaled  to 
the  slippage  distance  A/A.,  the  distance  the  electron  pulse  slips  back  relative  to  the 
optical  pulse  over  one  pass  through  the  undulator.  This  is  the  distance  over  which  the 
optical  field  and  the  electron  beam  interact.  All  longitudinal  distances  are  normalized 
to  this  slippage  distance  so  that  z/A/A  z.  This  then  gives 

***  FEL  Phase  Space  Evolution  *** 
j=932  aQ=0  Vq=0  N=413 
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Figure  3-7.  Phase-space  plot  for  40  A  X  ray  FEL, 
(With  dispersive  section). 
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a  a(z),  v  v(z),  and  C  C(2)  over  a  number  of  sites.  The  pendulum  and  wave 
equation  can  be  generalized  to  [3] 


^^z|COS(^^^  +  (j)^) 


(3.1) 


and 


^^2  •  —i  t 

-^  =  -yz«<®  .  (3.2) 

where  the  longitudinal  field  sites  z  refer  to  a  position  in  the  optical  field  envelope. 
When  the  electrons  pass  through  the  undulator  they  exchange  information  with  the 
optical  wave  envelope  in  a  range  of  sites.  The  optical  light  traveling  at  speed  c. 
remains  fixed  in  z.  The  slower  electrons  slip  back  to  site  z  - 1.  Those  electrons  at 
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site  z  interact  with  a  range  of  sites  in  the  optical  wave  envelope.  The  electrons  then 
pass  information  from  one  site  to  another. 

Longitudinal  multimode  simulations  shown  in  Fig.  3-8  and  3-9  using  the  same 
parameters  as  in  Fig.  3.7  and  3.8,  study  the  coherence  development  of  the  optical 
field.  At  the  bottom-left  of  each  simulation  we  see  the  gain  over  the  undulator  length. 
The  bottom-middle  is  the  single-mode  gain  spectrum  as  a  function  of  frequency 
centered  around  resonance.  The  bottom  right  is  the  evolution  of  power  developed 
through  the  undulator  length,  measured  in  a  dimensionless  power,  corresponding  to 
tens  of  GW.  The  middle-left  plot  is  the  growth  in  the  optical  field  amplitude  over  two 
slippage  distances  in  length,  The  middle  and  top-middle  show  the 

evolution  of  the  power  spectrum.  The  triangle  tick  mark  indicates  resonance,  v  =  0 
and  the  rectangle  tickmark  indicates  v  of  the  electron  beam  at  the  end  of  the 
undulator.  In  the  klystron  undulator,  the  final  power  is  much  higher.  The  middle-right 
and  top-right  frames  show  the  development  of  the  electron  phase  velocity  distribution. 
There  is  a  decrease  in  the  average  electron  phase  velocity  which  corresponds  to 
energy  given  up  to  the  optical  field.  In  the  FEL  without  a  klystron,  little  energy  is  given 
up  to  the  optical  beam,  while  the  system  with  a  klystron  gives  up  significantly  more 
energy.  The  simulation  without  a  klystron  and  before  saturation  resulted  in  an  X-ray 
linewidth  of  AA/X  «  0.1%,  where  AA/A  =  Av/27iA/  and  H  is  the  number  of  periods  in  the 
undulator.  Adding  a  dispersive  section  does  not  change  the  linewidth.  There  is  no 
degradation  to  the  quality  of  the  optical  beam. 

In  these  systems  the  standard  deviation  of  phase  velocity  was  oq  -  1 .04.  As  the 
standard  deviation  of  the  phase  velocity  inaeases  it  makes  it  difficult  for  the  electrons 
to  bunch.  Tin  phase  velocity,  oq,  must  be  considered  when  using  a  klystron  FEL. 
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******  FEL  Wrap  Evolution  ****** 
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Figure  3-8.  Longitudinal  multi-mode  simulation  for  40  A  X  ray  FEL, 
(Without  dispersive  section). 
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Rgure  3-9.  Longitudinal  multi-mode  simulation  for  40  A  X  ray  FEL, 
(With  dispersive  section). 
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C.  PHASE  VELOCITY  SPREADS  AND  KLYSTRONS 


As  mentioned  in  Section  II.E,  for  a  change  in  initial  energy  Aymc^,  there  is  a 
corresponding  change  in  electron  phase  velocity.  Av  »  AnNAyly,  in  which  the  Gaussian 
distribution  of  energies  with  an  rms  spread  .  Ayly,  characterizes  the  standard  deviation 
oq  =  AnNAyly.  When  oq  =  n,  bunching  becomes  impaired.  The  use  of  a  klystron  to 
bunch  the  electrons  may  not  have  a  positive  effect  if  the  electron  energy  spread  is  so 
large  that  the  modulator  cannot  develop  the  phase  velocity  differences  to  prepare  the 
electrons  for  bunching.  Figures  3-1  Oa.  3-1  Ob,  3-1  Oc.  and  3-1  Od  are  plots  on  how  the 
standard  deviation  in  energy  spread  can  effect  the  use  of  a  klystron  undulator  at  three 
different  dimensionless  current. 

Figures  3-1  Oa,  3-1  Ob.  and  3-1  Oc  show  what  happens  with  various  energy 
spreads  on  high-current  systems  {j  =  10000,  j  =  2200,  and  y  =  500).  Rgure  3-1  Od  is 
the  plot  for  low-current  system  (y  =  1)- 

In  all  four  figures  you  can  see  as  the  energy  spread,  oq,  inaeases,  the  overall 
gain  of  the  system  decreases.  By  using  a  klystron  undulator  it  is  hoped  to  increase 
the  gain  by  at  least  double  the  magnitude  compared  to  the  case  where  a  klystron  was 
not  used.  With  a  high  oq.  it  is  apparent  from  the  figures  that  a  klystron  will  not 
produce  a  gain  significant  enough  to  warrant  its  use.  The  use  of  a  klystron  at  the  high 
Oq  (2.0  and  3.14)  may  even  cause  a  decrease  in  gain.  In  the  monoenergetic  beam 
(ctq  =  0.0),  there  is  no  limit  to  the  strength  of  the  dispersive  section  or  to  the  gain  of 
the  system;  but,  this  is  unrealistic.  The  plots  indicate  that  a  moderate  spread, 
0.0  <  Oq  <  1 .5.  can  be  compensated  for  in  a  klystron  and  provide  a  significant  gain 
that  would  warrant  its  use.  The  magnitude  of  D  is  critical  though.  If  D  is  to  large  it 
may  have  a  negative  effect  and  cause  the  system  to  decrease  in  gain.  It  is  also 
shown  that  as  the  dimensionless  current  is  increased  the  critical  size  of  D  also 
increases.  A  lower  y  can  mean  a  lower  D,  the  distance  for  the  electrons  to  travel  to 
bunch  adequately  is  smaller  in  the  lower  energetic  beam. 

It  was  noted  that  for  the  4.5  A  and  the  1 .5  A  lasers  that  SLAG  has  proposed  the 
initial  electron  phase  vetodty  has  a  standard  deviation  of  oq  =  2.2  and  4.2.  The 
strength  of  the  dispersive  section  thus  becomes  critical. 
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D.  1 .5  A  FEL  WITH  A  KLYSTRON 


The  use  of  a  dispersive  section  to  shorten  the  length  of  a  high  power  FEL  such 
as  proposed  at  SLAG  would  reduce  the  cost  of  the  system.  The  linac  is  already  in 
place,  therefore  the  cost  to  build  the  FEL  is  largely  in  the  manufacturing  of  an 
undulator.  Shortening  the  undulator  would  reduce  its  cost  significantly  since  the 
manufacturing  ofundulator  sections  is  large. 

To  achieve  significant  gain  in  the  1 .5  A  FEL.  the  strength  of  a  dispersive  section 
for  y  =  3121  and  oq  =  4.2  is  estimated  to  be  0.2  from  Rg.  3-1  Ob.  With  high  current, 
the  maximum  gain  of  an  FEL  is  obtained  when  vq  is  slightly  off  resonance.  When 
adding  a  dispersive  section  to  the  undulator,  the  maximum  gain  is  shifted  to  a  higher 
value  of  Vq,  where  vq  =  nl2D.  The  optimum  for  this  system  is  found  to  be  vq  »  4.0. 

Figure  3-1 1  is  the  phase-space  plot  for  a  1 .5  A  FEL  with  a  length  of  30  m,  a 
dispersive  strength  of  D  =  0.2,  j  -  2170,  gq  =  3.8,  and  vq  =  4.  The  length  of  the  FEL 
was  shortened  from  35  m  to  30  m  which  in  turn  lowered  j  and  gq.  The  phase-space 
plot  shows  the  system  saturating  at  the  end  of  the  undulator.  The  peak  power  out  is 
14  GW.  which  is  slightly  higher  than  the  simulation  run  without  a  dispersive  section. 
All  indications  show  that  a  dispersive  section  can  be  effective  in  shortening  the  length 
of  the  undulator  and  increasing  gain. 

In  order  to  examine  the  coherence  development,  we  look  at  a  longitudinal  multi- 
mode  simulation  in  Rg.  3-12.  We  can  see  in  the  middle  window  a  narrow  spike 
indicating  a  narrow  band  width  of  AX/X  =  .1%,  which  indicates  good  coherence 
development.  The  right  window  shows  an  average  decrease  in  electron  phase  velocity 
which  corresponds  to  the  electrons  giving  up  energy  to  the  optical  beam,  creating  high 
optical  power. 

All  the  simulations  that  were  run  show  that  a  dispersive  section  may  have  a 
positive  affect  in  shortening  the  length  of  the  undulator.  The  use  of  a  dispersive 
section  may  decrease  the  cost  of  the  FEL  to  be  built  by  SLAG. 
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Rgure  3-1 1 .  Phase-space  plot  for  1 .5  A  system,  undulator  length  of  30  m. 

(With  dispersive  section) 


Figure  3-12.  Longitudinal  multi-mode  simulation  for  1.5  A  X  ray  PEL. 
(With  dispersive  section). 
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IV.  CONCLUSIONS 


The  theory  that  has  been  presented  is  not  limited  to  X  ray  FELs,  but  to  all  FELs. 
The  Free  Electron  Laser  has  many  uses  in  the  military  and  industrial  applications: 
ranging  from  high-power,  speed-of-light  self-defense  weapon  systems  to  DNA-pair 
imaging.  The  ability  of  an  FEL  to  be  tunable  allows  it  to  overcome  limitations  that 
other  high  powered  lasers  face,  such  as  atmospheric  effects.  The  FEL  as  a  self 
defense  weapon  system  does  not  require  replenishment  from  a  limited  magazine.  The 
FEL  is  a  system  that  is  needed  in  the  present  day  military  to  overcome  conditions  that 
have  outgrown  systems  now  in  use.  There  is  a  need  for  the  United  States  to  be  the 
leader  in  laser  technology  to  ensure  its  capability  to  improve  world  conditions. 

The  proposed  SLAG  X-ray  FEL,  which  will  come  on-line  in  the  near  future,  will  be 
the  first  coherent  X  ray  light  source.  The  acquired  technology  is  not  limited  to  the 
SLAG  proposal,  but  can  be  used  to  develop  other  types  of  high  powered  FELs.  An 
injector  system  now  under  development  for  the  SLAG  linac  will  enable  it  to  produce  a 
4.5  A  or  1 .5  A  FEL. 

The  use  of  a  dispersive  section  will  aid  in  shortening  the  length  of  the  undulator. 
A  significant  part  of  the  cost  to  build  the  FEL  at  SLAG  is  manufacturing  the  undulator, 
so  shortening  of  the  undulator  is  economically  important.  However,  beam  quality  must 
be  high  in  order  to  use  a  dispersive  section.  So  the  cost  of  improving  the  injected 
electron  beam  quality  must  be  considered  in  the  decision  to  use  a  klystron  undulator. 
The  better  the  quality  of  the  electron  beam  the  greater  the  positive  affects  a  dispersive 
section  may  have  on  increasing  gain  and  shortening  the  overall  length  of  the  system. 

Finally,  we  have  presented  simulation  results  to  show  that  use  of  dispersive 
sections  is  a  possible  and  potential  cost-saving  option.  Further  investigation  in  the  use 
of  dispersive  sections  on  high-gain  FELs  is  required.  For  example,  positioning  the 
dispersive  section  along  the  undulator  to  optimize  its  effects  may  be  studied. 
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